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Abstract: We present numerical simulation of bulk Lorentz factor of relativistic electron- 
positron jet driven by Compton rocket effect from accretion disk radiation. The plasma 
is assumed to have a power-law distribution n e ( 7 ) oc 7 _s ( 1 < 7 < 7 max ) and is continu¬ 
ously reheated to compensate radiation losses. We include Klein-Nishina (hereafter KN) 
corrections, and study the role of energy upper cut-off 7 max, spectral index s, and source 
compactness. We determine terminal bulk Lorentz factor in the case of supermassive 
black holes relevant to AGN and stellar black holes relevant to galactic microquasars. In 
the latter case, effects of KN corrections are more important, and induce terminal bulk 
Lorentz factor 7 ^ smaller than in the former case. The result can explain the low bulk 
Lorentz factors for galactic sources (GRS1915+105, GROJ1655-40) compared to extra- 
galactic ones. We also take into account the influence of the size of the accretion disk; if 
the external radius is small enough, the bulk Lorentz factor can be as high as 60, which 
is comparable to the values needed to explain extragalactic gamma-ray bursts. 


1 Introduction 

Superluminal motion observed in active galactic nuclei (AGN), especially in the 
class of blazars, seems to be closely linked with high-energy emission. Such motion 
was recently observed in the Galaxy (Mirabel & Rodriguez, 1994, Hjellming & 
Rupen, 1995, Tingay, S. J., et a/., 1995) in the so-called microquasars. Nevertheless 
differences are underlined in those two cases. The latter systems were observed with 
small value of bulk Lorentz factor (around 2.5), while in the former ones values 
of about 10-20 are frequent. A possible acceleration mechanism is the so-called 
’Compton rocket’ effect (O’Dell, 1981), i.e. anisotropic inverse Compton effect on 
a highly relativistic plasma. R was then argued (Phinney , 1982) that because of 
Compton cooling, only small value of 7^00 could be reached by this mechanism. 
However, in the frame of the ’two-flow’ model (Sol, Pelletier & Asseo, 1989), a pair 
plasma can be reheated by the turbulence triggered by a surrounding jet (Henri & 
Pelletier, 1991) and then the Compton rocket can be more efficient (Marcowith et 
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al ., 1995). We here describe more precisely our assumptions and present our main 
results on acceleration of such pair plasma. 


2 Compton rocket effect with Klein-Nishina corrections 


The pair plasma is assumed to be described by an energy distribution n' e (z, Y) oc 
7 /-s for 7 min < l' < Imax , with s, 7 m in and -frnax independant on 2 . The radiation 
force is due to soft photon coming from a standard accretion disk (Shakura & 
Sunyaev, 1973) around a Schwarzschild black hole. 

Let be the rate of momentum change for one electron and /3'c— yj 1 —i I/ 7 /2 c 
its velocity. The force exerced upon the pair plasma is then in the blob rest frame: 


F' z = j df2' e d'y , n' e ('y', f2' e )^-.u z , ( 1 ) 

where u z defines the jet axis direction. To determine F ,z . one can use two 
different approximations: for low energy particules, the KN cross-section tends to 
the Thomson limit and we can apply the result of Phinney (1982) for an isotropic 
pair plasma. For high energy particules, KN corrections become important but 
we can assume that the head-on approximation is valid. In this approximation we 
determine both momentum and energy loss rates in the electron rest frame and rely 
them to the blob rest frame using Lorentz transformation (see Renaud & Henri, 
1997, for detailed calculations.). We connect the two regimes by defining a critical 
Lorentz factor Y crit for which errors in the head-on approximation ( 7 ' > 7 cr it) and 
in the Thomson limit ( 7 ' < 7 cr it) are of same order. This gives 7 cr it ( 2 (e))- 1 / 3 . 
(e) oc A'f 1 / 4 M -1 / 2 , is the average photon energy emitted from accretion disk, M is 
the accretion rate and M the mass of the central black hole. For AGN (s) ~ 10 -4 
and q /era ~ 20 — 30 while for a microquasar (e) ~ 10 -2 and 7 cr i t ~ 5. Following 
Phinney (1982) we determine the acceleration of pair plasma by considering the 
conservation of stress-energy tensor. For a reheated relativistic plasma (with p r — 
p'/ 3) one finds the equation of motion. 


d^b F' z 1 

dz p' Aj- + 1 ’ 
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3 Results 

3.1 Effect of Klein-Nishina corrections 

The effect of the KN corrections is to reduce the contribution of high energetic 
collisions in the net radiation force seen by the plasma. This leads to smaller value 
of terminal Lorentz factor than expected in the Thomson regime. KN corrections 
begin to dominate roughly when 7 max{F) ~ 1- The importance of this effect is 
so predominant for microquasars in which accretion disk radiates more energetic 
photons than in the extragalactic case. 
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3.2 Influence of the different parameters 


We determine terminal Lorentz factor for different configurations: we consider the 
case of supermassive black holes M — 10 9 Mq relevant to AGN and stellar black 
holes M — 5Mq relevant to microquasars. For these two typical cases we show 
figure 1 the influence of the spectral index and the energy cut-off 7 max . 
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Fig. 1. Terminal Lorentz factor 'fbca as a function of y maa; for different value of spectral 
index s. The two top pannels correspond to M — 10 9 Mq black hole with r e = 3.10 3 r ff , 
L = Lfidd (left) and L = O.lLEdd (right). The two bottom pannels correspond to M = 
5 Mq black hole with r e — 3.10‘Vg, L — L E dd (left) and L — 0.1L E dd (right). 


For s < 3, F ,z is dominated by the high value of y 7 . Terminal Lorentz factor 
admits then a maximum when KN corrections begin to dominate. This value is 
more rapidily reached for a stellar black hole, explaining smaller 7 ^ than in 
extragalactic configurations. For s > 3, F ,z is dominated by the lower cut-off of 
the electronic distribution, explaining why no variation is evident with 7 maa; in this 
case (see Renaud & Henri, 1997). The plasma behaves as a cold one, and one finds 
small value of terminal Lorentz factor, like Phinney’s result. The dependance on r e 
is shown figure 2, where we extend the previous calculation to the case r e = lOry;. 
We find the same global behaviour of 'fboo with spectral index and 7 ma x- In this 
configuration the radiation force is more efficient and can be as high as 60 in 
the extragalactic case. 






Fig. 2. Terminal Lorentz factor as a function of 7 max f° r different value of spectral 
index s. The two pannels correspond to r e = 10 r g and L = LEdd, M — 10 9 M©(left) 
M — 5 Mq (right). 


4 Conclusion 


We find that for a given luminosity, the terminal Lorentz factor 'ytoo admits a 
maximum due to KN corrections. They are predominant for stellar black holes 
leading to smaller values than for a supermassive black hole. Moreover the observed 
value of Lorentz factor of about 2.5 for the two microquasars with large spectral 
indexes (respectively s ~ 4, Finoguenov et al, 1994 and s ~ 4.6, Harmon et al, 
1995) are coherent with our results. We also find that for fixed mass, spectral index, 
7 boo increases with decreasing r e , reaching 60 for r e = 10r, y , which is comparable 
to the value needed to explain observed flux in gamma-ray burst. 
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